Parametric Resonance of Optically Trapped Aerosols 
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The Brownian dynamics of an optically trapped water droplet are investigated across the transi- 
tion from over to under-damped oscillations. The spectrum of position fluctuations evolves from a 
Lorentzian shape typical of over-damped systems (beads in liquid solvents), to a damped harmonic 
oscillator spectrum showing a resonance peak. In this later under-damped regime, we excite para- 
metric resonance by periodically modulating the trapping power at twice the resonant frequency. 
The power spectra of position fluctuations are in excellent agreement with the obtained analytical 
solutions of a parametrically modulated Langevin equation. 



Parametric resonance provides an efficient and 
straightforward way to pump energy into an under- 
damped harmonic oscillator In general, if the res- 
onance frequency of an oscillator is dependent upon a 
number of parameters modulating any of these at twice 
the natural oscillation frequency parametrically excites 
the resonance. Such behavior leads to surprising phe- 
nomena in the macroscopic world (pumping a swing, sta- 
bility of vessels, surface waves in vibrated fluids) [1, 
On the microscopic scale, where stochastic forces become 
important, one refers to Brownian parametric oscillators 
0]. As an example, the parametric driving of Brownian 
systems has been shown to be at the origin of some pe- 
culiar behaviors such as the squeezing of thermal noise 
in Paul traps [5[. Parametrically excited torsional oscil- 
lations have also been reported in a single-crystal silicon 
microelectromechanical system What makes para- 
metric resonance useful is that in many cases it is easier 
to modulate a system parameter rather than applying 
an oscillating driving force. Moreover, for finite but low 
damping rates, we may never reach a stationary state 
with the damping forces dissipating all of the input power 
and consequently the amplitude of oscillations diverge. 
Optically trapped microparticles constitute a beautiful 
example of Brownian damped harmonic oscillator (DHO) 
and they are becoming an increasingly common tool for 
the investigation of different fields of basic and applied 
science . The possibility of pumping mechanical energy 
into optically trapped particles could open the way to 
many applications. In optical tweezers, even though it is 
easy to periodically modulate the laser power, paramet- 
ric excitation is usually ineffective because of the heavy 
damping action of the surrounding fluid. 

Recently it has been reported that modulating the 
laser power at the parametric resonant frequency in an 



overdamped system increased the amplitude of mean 
squared fluctuations [8|. However, these findings have 
been difficult to reproduce and are in strong contrast 
with the prediction of Langevin dynamics [1, [TO, [lH • 

In this Letter we show how parametric resonance can 
be excited in optically trapped water droplets suspended 
in air, due to the reduced damping force. We measure 
power spectra of position fluctuations and find an excel- 
lent agreement with the theoretical expectations based 
on Langevin dynamics with a parametric forcing. 

The dynamics of an optically trapped droplet is de- 
scribed by the Langevin equation 
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where = k/m is the natural angular frequency of 
the oscillator depending on trap stiffness k and particle 
mass m. Tq = 67Tr]a/CcTn is the viscous damping due to 
the medium viscosity r] and depending on particle radius 
a and mass m. One must be careful to consider the non- 
continuum effects of Stokes' Law in air due to the finite 
Knudsen number of the particles under study. To correct 
Stokes' Law the empirical slip correction factor, Cc, is 
introduced, with a 5.5-1.6% reduction in drag for 3-10/im 
diameter droplets respectively [I^. 

The stochastic force ^, due to thermal agitation of sol- 
vent molecules, is generally assumed to be uncorr elated 
on the time-scale of particle's motion: 



(mm) = 2ToKBT/m 5{t) 



(2) 



The corresponding power spectrum S^{(jo) of the 
stochastic variable ^ can be defined as: 
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where: 



i{uj) = l/2TT ^{t) exp[-iujt]dt 
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By putting (|4j) and ([2]) in (|3]) it is easily shown that ^ 
has a white noise spectrum 



TqKbT 



(5) 



By Fourier transforming ([T]) and using ([5]) we can easily 
obtain the power spectrum of position fluctuations as: 
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When trapping objects in liquid solvents the ratio 
Vq/VLq is always larger than one, i.e the system is over- 
damped. The ratio Tq/VLq depends only slightly on par- 
ticle's size and, in solvents with water-like viscosities, is 
always greater than 1 up to power levels of some tens 
of Ws. For typical trapping powers of order 10 mW in 
water Fq/I^o is typically > 10. As a result only those fre- 
quencies smaller than Fq, and hence much smaller than 
r^o, have a significant amplitude in the power spectrum 
of X. Under these conditions we can therefore neglect cj^ 
with respect to Vt^ in the first term of the denominator 
in (|6|) and obtain the usual Forentzian power spectrum 
characterised by o;"^ tails [13]. Such an overdamped 
condition precludes the possibility of exciting significant 
oscillations either directly or parametrically. To probe 
oscillations in the liquid damped regime we would need 
to be able to increase typical trap power by four orders 
of magnitude - introducing uncontrollable heating and 
damage of the trapped object. A more feasible route is 
to reduce viscosity by two orders of magnitude. This last 
condition can be readily obtained by trapping particles 
in air whose viscosity is approximately 1/55^^ of water 
[t] = 1.8 X 10-^ Pa s) [13. 

For these experiments, our optical tweezers are based 
around an inverted microscope with a high numerical 
aperture oil immersion microscope objective (1.3NA, 
100 x). The continuous wave laser is a Nd:YAG, fre- 
quency doubled to give 0-2 W of 532 nm light. To couple 
the beam into the air medium, a single cover slip is rested 
over the objective on a thin oil layer. A water aerosol is 
produced using a nebulizer and usually a 3-10/im diam- 
eter water droplet is trapped at the beam focus lH, 
see Fig. 1. One should note that whereas for particles 
trapped in fluid, a laser power of 10s of mW is typical, 
here, to maximise the stiffness of the trap, we use powers 
of order 1 W. 



FIG. 1: The measured power spectra of trapped aerosol par- 
ticle at two different powers. At lower power (black circles) it 
is overdamped and the mean squared amplitude of the high 
frequency motion decays as cj~^. At higher powers (white 
circles), the aerosol is underdamped and the mean squared 
amplitude decays as o;"^. The inset shows an optical image 
of a trapped aerosol particle. 



A quadrant photo detector, placed in the back focal 
plane of the condenser lens, receives the light transmitted 
through the droplet. By measuring the imbalance of the 
light collected by the quadrants, the lateral displacement 
of the droplet is deduced with a bandwidth of several kHz 
and a precision of better than 5 nm [l7i]. 

The power spectra of the measured displacement, for 
two different trap powers, are shown in Fig. [TJ It is 
clearly visible how the particle dynamics crosses over 
from an over damped dynamics with a Forentzian spec- 
trum with a high frequency roll-off proportional to o;"^ 
to an underdamped regime with a faster roll-off, uo~^^ 
and the appearance of a resonance peak at a frequency 
of about 1 kHz. The emergence of such a peak arises from 
the fact that the inertial terms in ([1]) are no longer neg- 
ligible. As a consequence an average trajectory starting 
away from the equilibrium position crosses the equilib- 
rium position with a finite velocity. 

In this situation the parametric resonance is excited 
by modulating the strength of the trapping potential. 
Ideally the potential is made shallower when the particle 
traverses the equilibrium position and steeper again when 
the particle is far from the equilibrium position. This is 
maximally efficient when we modulate the potential at 
twice the natural oscillation frequency Vt^. To consider 
this model in quantitative terms we can rewrite ([1]) in 
the presence of a parametrically modulated external po- 
tential: 
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trap. The white circles of Fig|2]show the measured power 
spectrum for a trapped water droplet when the trap has 
constant power. The presence of the peak suggests that 
we are in an under damped regime. By fitting to equation 
(|6]) (solid line) we can directly extract the resonant fre- 
quency l]o/27r = 2.0 kHz and the damping term To = 6.8 
kHz relevant to our experimental conditions. The fitted 
value of To corresponds to the Stokes drag on an aerosol 
droplet of radius 3.4/im. 

We then apply the square- wave modulation of the trap- 
ping power, with the laser adjusted to give the same av- 
erage power as before, fti :^ 21^0 and g = 0.4. The black 
circles of Fig|2] show the power spectrum of the lateral 
motion for the trapped water droplet with modulated 
laser power. The excitation of a resonance appears as a 
higher and narrower peak in the power spectrum. This 
observed spectrum matches closely the expected behavior 
(solid line) obtained by applying the measured parame- 
ters r^o, To, r^i, g to equation ([T4j), strongly supporting 
our interpretation of this peak being due to a parametric 
excitation of the resonance. Higher harmonics, character- 
izing the response of parametrically driven systems, are 
suppressed in our case, being only slightly underdamped. 

Using these parameters with (p!4|) we can make gen- 
eral predictions about the dynamics of a particle that 
can then be verified against experiments. One compar- 
ison to make is the predicted and observed form of the 
power spectra as a function of the modulation frequency, 
both above and below the parametric resonance condi- 
tion. These results are shown in Fig. [3l the white circles 
are the experimental points and the black lines are the 
predicted spectra. Again, there is an excellent agreement 
between the observed and predicted particle motion. In 
particular the parametric excitation of oscillations man- 
ifests as a narrowing of the peak (or a reduced apparent 
damping F, defined as the full width half maximum), oc- 
curring when modulating at twice Qq. A shift in peak po- 
^nki^) = [-{^^nnif^nl^i{u;^nni)To]Snk^^lgak-n ^^^^^^ ^^^^ apparent close to parametric resonance 



x{t) -^nlil^ gf{t)]x{t) + Fox(t) = at) (7) 
/(t + T)=/(t), -l<f{t)<l (8) 

where < ^ < 1 measures the strength of modulation. 
By Fourier transforming ([7j) we obtain: 



{—ij'^-\-Ql-\-iLjro)x{Lj)-\-QQg ^ akx{ij -\-kQi) = ^{ij) 

k= — oo 

(9) 

where is the coefficient of the /c27r/T = kQi fre- 
quency component of the Fourier series expansion of f{t). 
It is clear from equation (|9|) how parametric modulation 
introduces a coupling between all those frequencies dif- 
fering by an integer number of Qi. We now introduce the 
vectors Xn{uj) = x{uj -\- nQi) and Rn{^) = + nQi) 
and write the recursive relations: 



[-{uj + nnif ^nl^ i{uj + nni)ro]Xn{uj) + 

oo 

^ akXn^ki^) = Rfi' (10) 

k— — oo 

To obtain the power spectrum Sx{oo)^ for each fre- 
quency UJ we should compute Xo{uj). This will be, in 
turn, coupled to all other components in the array X^. 
However the strength of the coupling will decay for large 
\n\ so that we can limit ourselves to a finite number of 
components and write the matrix equation for the array 



G-i(a;)XH =R(a;) 
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with 



(12) 

By matrix inversion we obtain the power spectrum as: 
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and from the definition of power spectrum: 
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If fto and Fq are known, we can use to predict the 
power spectrum of a Brownian particle in a modulated 



Both of these signatures are compared with theoretical 
expectations in Fig. [4] further supporting our interpreta- 
tion of the system as being a Brownian parametric oscil- 
lator . 

We recognise that the system we report here relates to 
study of only the lateral motion of the trapped droplets. 
In keeping with other work we note from examination 
of the video images that the axial movement of the 
trapped droplet has significant amplitude on much longer 
timescales, corresponding to frequencies in the region of 
10-50Hz [19|. This reflects the comparatively weak ax- 
ial trapping, possibly arising from aberrations associated 
with non optimised objectives or the increased scattering 
force. It may be possible to use doughnut or Laguerre- 
Gaussian beams having zero on-axis intensity, and im- 
proved axial trapping [l8[ . 

We have reported the first observation of a paramet- 
rically excited resonance within a Brownian oscillator. 
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FIG. 2: The measure power spectrum of a trapped water 
droplet for no modulation of the laser power (white circles) 
and modulation at 3.9 kHz (Qi ^ 2Qo) (black circles). The 
peak is higher and narrower on the resonant condition thus 
indicating parametric excitation. Solid line below black circles 
is the predicted spectrum from ([14]). 
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FIG. 3: Evolution of position power spectra on varying 
the modulation frequency VL\. Parametric excitation of os- 
cillations is evident at the parametric resonance condition 
^i/^o = 2. The solid lines are the theoretical predictions 
from (dH). 



The demonstration of this effect within optical tweezers 
for this purpose was made possible by relying on the vis- 
cosity of air to lightly damp the motion of a trapped 
aerosol droplet. The detailed observed dynamics match 
closely the power spectra predicted from a parametrically 
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FIG. 4: a) Peak fullwidth half maximum as a function of 
modulation frequency. Solid lines are the theoretical predic- 
tions from ()14|) . b) Peak position shift as a function of mod- 
ulation frequency. 



modulated Langevin equation, where all parameters were 
fixed from examination of the non modulated system. 
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